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Abstract 



We investigate the prospects of atomic interference using samples of Bose 
condensed atoms. First we show the ability of two independent Bose conden- 
sates to create an interference pattern. This holds even if both condensates 
are described by Fock states. Thus, the existence of an experimental signa- 
ture for a broken gauge symmetry, seen in a single run of the experiment, is 
not necessarily reflected by a broken symmetry on the level of the quantum 
mechanical state vector. Based on these results, we simulate numerically a 
recent experiment with two independent Bose condensates [K.B. Davis et al., 
PRL 75, 3969 (1995)]. The existence of interference fringes is predicted based 
on the nonlinear Schrodinger equation. Finally we study theoretically the in- 
fluence of finite temperatures on the visibility of the interference in a double 
pinhole configuration. 

03.75.Fi,05.30.Jp 
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I. INTRODUCTION 



Recent realizations of Bose-Einstein condensation (BEC) in dilute and ultracold gases 
]l]j2||J have attracted vivid interest. It is hoped that the study of those experimental systems 
will give new insight into the physics of BEC. Since the current understanding of BEC is 
largely influenced by the concept of a macroscopic wave function, the study of this feature 
is of foremost importance. The investigation of interference phenomena should be perfectly 
apted for this purpose. Another motivation for the study of interference properties is the en- 
visioned development of a new source of atoms, based on BEC, with high flux and coherence, 
that is expected to stimulate atomic interference experiments. 

Although interference is a well known phenomenon in quantum statistically degenerate 
systems (cf. Josephson effect), surprisingly little quantitative analysis has been performed 
up to now, to investigate the interference properties of a Bose condensate. In this paper we 
study theoretically two interference experiments that could be performed with the technology 
of current BEC experiments 

In Sec. [H]we focus on the general possibility of interference between two independent Bose 
condensates. This topic is closely related to earlier studies, showing interference between 
two independent laser beams Interference of atoms, originating from independent 

sources, differs qualitatively from usual interference experiments, where atoms interfere only 
with themselves. In this case the phase of the interference pattern is uniquely defined by 
the geometry of the setup. For two independent condensates, however, an interference 
pattern will be exhibited with a phase that unpredictably varies between different runs of 
the experiment. This implies that the system possesses some symmetry that is spontaneously 
broken in a single run of the experiment. An intuitive explanation for this interference follows 
from the common notion of spontaneously broken gauge symmetry, implying the existence of 
a macroscopic wave function. We show, however, that the system behaves identically under 
circumstances where the broken symmetry is not reflected by the quantum mechanical state 
vector, e.g. if the two condensates are initially in Fock states. In this case, it is still possible 
to derive from quantum mechanics the existence of interference patterns in a single run. 

Having shown the interference of two independent Bose condensates in quite general 
terms, we turn more realistic in Sec. |T|. This includes accounting for the finite size of exper- 
imental condensates and the important influence of atom interactions. Finite temperature 
effects, however, are still neglected. In principle, a setup for interference of two condensates 
has been realized in the BEC experiment of || . There, two condensates are stored in a mag- 
netic trap, separated by a light beam. After releasing them from the trap, they expand and 
eventually overlap. However, no interference has been observed in the reported experiment 
||. We simulate this experiment numerically using the nonlinear Schrodinger equation. We 
are able to reproduce the observed final size of the atomic cloud in good agreement with || . 
In addition it is shown that interference fringes should exist with a spatial period smaller 
than the present experimental resolution. 

In Sec. [TV] we investigate the more conventional case where a single condensate is released 
from a harmonic trap and subsequently interferes with itself in a simple double pinhole ex- 
periment. In this calculation we allow for finite temperatures but neglect atom interactions. 
It is investigated, to what extent the existence of interference fringes might serve as a sig- 
nature of Bose condensation. It turns out that the appearance of a condensate is reflected 
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by an extremely sharp change in the visibility of the interference pattern at the critical 
temperature. 

After completion of our work we have become aware of a recent paper by J. Javanainen 
and S. Yoo ||. These authors have studied the interference of two independent Bose conden- 
sates with conclusions, similar to the ones reached in Sec. [II], but based on a complementary 
derivation. 

II. INTERFERENCE OF TWO INDEPENDENT BOSE CONDENSATES 
A. Spontaneously broken gauge symmetry 

In usual interference experiments (cf. Sec. [TV]) , atomic beams are split apart in a suitable 
way, e.g. by a Young double slit, and recombined afterwards thereby leading to an interfer- 
ence signal. Since all atoms come from the same source and are not tracked on their way 
through the interferometer one usually notes that atoms interfere with themselves in such a 
situation. Here, we address the question whether an interference signal can be exhibited as 
well, when two independent Bose condensates merge. 

Although, in this case atoms do not interfere with themselves but rather with other 
indistinguishable atoms, it should indeed be possible to observe such an interference under 
suitable conditions. A simple understanding of this phenomenon can be achieved using the 
notion of spontaneously broken gauge symmetry, that is widely regarded as a characteristic 
feature of BEC. We will shortly summarize this concept in the following. 

A spontaneously broken symmetry, in general, implies that the behaviour of a single 
many-particle system differs from its ensemble average. The common paradigm for sponta- 
neous symmetry breaking is a ferromagnet, where a single domain may show a magnetization 
while the mean magnetization, averaged over different domains, vanishes. It is a common 
approach in statistical physics [0] to investigate the existence of such a broken symmetry by 
adding a suitable small symmetry breaking field (Bogoliubov auxiliary field) to the Hamil- 
tonian of the system. The ground state of this new Hamiltonian is meant to describe a 
particular single system at zero temperature with broken symmetry instead of the symmet- 
ric ensemble. The density matrix, describing the full ensemble, is then assumed to be a 
symmetry preserving, incoherent superposition of degenerate ground states, corresponding 
to different orientations of the Bogoliubov field. 

So as to account for several observed interference phenomena like the Josephson effect, 
it is a common notion to postulate a spontaneously broken gauge symmetry || in the case 
of BEC. The corresponding Bogoliubov field consists in adding 

etf(x) +e*^(x) (2.1) 

to the Hamiltonian, letting e — ^ after having performed the thermodynamic limit. It 
does not correspond to any physical interaction, in contrast to the case of ferromagnetism. 
Its consequence is that a single Bose condensate is described by a coherent state with a 
distinguished phase. All measurable quantities of this particular single condensate are then 
fully characterized by the complex scalar field ip(x) with 

(*{x)) g = 1>{x). (2.2) 
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Note, that the angular brackets in Eq. ( |2.2|) refer only to the subensemble that is defined 
by the particular coherent state. The field ip(x) is usually called the macroscopic wave 
function of the condensate. The macroscopic wave function if)'(x) of any other condensate, 
seen in a different run of the same experimental setup, can be achieved by applying the 
gauge transformation 

¥(x) = exp(iaN)^(x) exp(-iaN) (2.3) 
with N being the number operator. This yields 

^'( x ) = e~ ia tlj(x), (2.4) 

motivating the notion of a broken gauge symmetry. In contrast to the single system, the 
full ensemble does not show any phase dependence 

(*(*)) = 0, (2.5) 

i.e. it is invariant under the symmetry transformation Eq. ( |2.3D . The main advantage of this 
concept is that quantities, measured in a single run of the experiment, can still be expressed 
by simple quantum mechanical expectation values, while the true ensemble expectation 
values fail to characterize a single system properly. 

So far, the analogy between spontaneous symmetry breaking in a ferromagnet and in BEC 
seems to hold perfectly. However, spontaneously broken gauge symmetry has to be taken 
with some caveat, as will be indicated in the following. In the case of ferromagnetism, the 
ground state of the unmodified Hamiltonian is degenerate. The effect of the Bogoliubov field 
consists only in selecting a particular one of these degenerate states for the characterization 
of a given domain. 

Such a correspondence exists equally in the usual treatment of BEC. Here, the density 
matrix of the condensate at T = is expressed in terms of coherent states 

p Q = - [ al 2 a 5{\a\ 2 - N) \a) (a\ , (2.6) 

7T J 

from which the Bogoliubov field selects a particular \a), due to the assumed broken gauge 
symmetry. 

The justification of this treatment may be questioned, since the Bogoliubov field Eq. 
( |2.1|) does not correspond to a real physical field. In addition, atomic coherent states are no 
eigenstates of the particle number and violate atom number conservation, that is implied by 
fundamental superselection rules. It is therefore not obvious why a decomposition in terms of 
coherent states rather than in terms of number states is appropriate for the characterization 
of a single run of the experiment. Another subtle feature arises from the fact that the 
assumed degeneracy of the ground state requires the use of the grandcanonical ensemble 
with fluctuations of the particle number. In a microcanonical ensemble with a fixed particle 
number N, however, the ground state of the unmodified Hamiltonian may be nondegenerate, 
i.e. given by a single Fock state 

Po = \N) (N\ . (2.7) 
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Obviously, the concept of the Bogoliubov field, merely selecting one of several existing ground 
states, becomes meaningless in this case. This raises the important question whether a 
microcanonical description will also show BEC with all its characteristic features. A positive 
answer to this question would imply that the degeneracy of the ground state is not essential 
for the understanding of BEC, in contrast to the case of ferromagnetism. Indeed, a numerical 
simulation of the quantum kinetic equation |J has shown that a macroscopic population of 
a nondegenerate ground state occurs in a system with a fixed particle number. It remains 
to be seen, whether the density matrix Eq. (2.7) is able to reproduce even those effects, that 
are commonly explained by a broken gauge symmetry. 

We address this question in the following by investigating the interference properties 
of two independent Bose condensates. It will be shown, that the notion of a macroscopic 
wavefunction gives the proper intuitive understanding about a single run of such a fictitious 
experiment, no matter whether Eq. ( [2.6|) or Eq. fl2.7|) is assumed. This indicates that a 
degeneracy of the ground state should not be essential for the explanation of the phenomena, 
that are usually regarded as an evidence for a spontaneously broken gauge symmetry. 



B. Interference of two plane waves 

We assume that two independent atomic beams, being Bose condensates, merge on a 
planar detection screen. The accumulation of a large number of particles on the screen 
during a given detection time will be regarded as a single run of the experiment. Only if 
we repeat the same experiment several times with a complete physical reset between different 
counting intervals, this will be called an ensemble. Since we assume a homogeneous velocity 
distribution of the incoming atoms, the distribution of particles, accumulated on the screen 
during in a single run is proportional to the line density N(x) along the x-axis. 

We describe the condensates by two independent macroscopically occupied momentum 
modes with x-momenta ±k and average occupation numbers iVj = (a\a^ (i = 1,2). Their 
number fluctuations are assumed to be negligible in comparison with the mean occupation 
numbers: 

alalia,) 1 + Q /ly 



N,N, \N, 

This includes a representation of the condensates in terms of Eq. (|2.6|) or alternatively of 

The first step of our analysis is to determine the ensemble behaviour of the system. It is 
easily seen from 

(N(x)) V 1/3 = N 1 + N 2 (2.9) 

that no interference is visible if the measurements of many runs are summed up. This result 
is no surprise due to the translational invariance of the system in x-direction. Usually, 
also a single run is sufficiently characterized by (N(x)). However, our considerations about 
broken gauge symmetry have already indicated that this may be wrong in the case of BEC. 
This implies that the density distribution N(x), measured in a single run, might show a 
nontrivial spatial dependence. The interesting question is therefore, whether interference 



5 



exists in the counts from a single run of the experiment. A rigorous answer can be given 
by analyzing correlation functions or conditional probabilites like (N(x)N(x + A)). If these 
show a nontrivial dependence on the relative distance A, it is clear that the individual 
density distribution N(x) must be different from the ensemble average Eq. ( |2.9|) . However, 
before turning to the detailed discussion of correlation functions, we describe briefly the 
predictions of the broken gauge symmetry model for our situation. 



1. Intuitive treatment 

Here, we assume the notion of a broken gauge symmetry to be true. As a consequence, 
the macroscopic wave function of a single condensate would be a classical plane wave with 
an arbitrary but fixed phase factor. If we describe the full state vector by a product of two 
coherent states, corresponding to the two condensates, we get the combined macroscopic 
wave function 

i){x)W T ^ = e iai ^N[e lh)X + e ia2 //V^T** *, (2.10) 

which is a sum of the single macroscopic wave functions. It is obvious that this wave function 
leads to an interference pattern 

N(x)V 1/3 = (Ni + N 2 ) + 2y / i\yV2 cos(2%r + <p) (2.11) 

on the screen which is determined by the relative phase ip = a% — a 2 - This relative phase 
varies between different runs and thus leads to Eq. ( |2.9| ) in the ensemble. The unpredictable 
variation of p between different runs would be usually interpreted as an experimental evi- 
dence for a broken gauge symmetry. 

However, the intuitive treatment using Eq. ( |2.10|) depends on assumptions, that were 
already questioned above. First, a degeneracy of the ground state is assumed, so as to de- 
compose the density matrix into coherent states (cf. Eq. ( |2.6|) ). Then the single condensates 
are identified with specific coherent states with the intention that their expectation value 
(ty^ (x)ty (x)j characterizes the density distribution of a single run. Regrettably, this pro- 
cedure is not in rigorous agreement with atom number conservation. These demand that 
the condensate of a single run has a definite atom number, although we may not know it. 
This number is allowed to fluctuate only between different runs of the experiment. Taking 
number conservation serious, one might be misled to suggest that the density distribution of 
a single run is properly characterized by the expectation value (ty^(x)ifl(x)^ of Fock states. 
This would yield 

N(x)V 1/3 = Ni+N 2 (2.12) 

even for a single run, which is in contradiction with Eq. ( |2.11| ). 

It will be the result of the rest of this section to confirm the implications of Eq. (|2.10|) 
rigorously without using any of the mentioned assumptions. A counterintuitive implication 
of this derivation is that Eq. Q2.ll ) instead of Eq. ( 2.12 ) gives the proper understanding of 



the single system behaviour even if the two condensates are described by Fock states. In 
addition, the following analysis illustrates, how one can extract detailed information about 
the behaviour of single systems from the quantum mechanical ensemble description. 
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2. Correlation function analysis 



Our aim, to characterize single runs of an experiment, confronts us with the problem that 
quantum mechanics by definition determines only properties of the ensemble, i.e. averages 
over many runs. Nevertheless, ensemble averages do allow to predict certain features of 
single runs which are not apparent in the ensemble by the analysis of correlation functions. 
Such an approach has found widespread use in quantum optics [[HJ. Characteristic features 
of the time evolution of a single atom could be derived in the case of the quantum jumps 



11,12], eventually leading to the numerical tool of the Monte-Carlo wavefunction technique 



13| , P^JT5| . Here, we consider the opposite limit where even a single run of the experiment 



involves a large number of particles. This implies that a complete spatial interference pattern 
can be built up in a single run. In our special situation we will be able to give a full 
analytical characterization of these interference patterns. An alternative approach would 
be to interprete the correlation functions in terms of conditional probabilities. Numerical 
simulations of single atom counts, based on these conditional probabilities, would lead to a 
similar result. 

In accordance with the experimental experience, we describe the density distribution 
N(x), measured in a single run of the experiment, as a smooth pattern N(x), blurred by 
some additional unpredictable and uncorrelated shot noise 

N(x) = N(x) +£(x). (2.13) 

The shot noise accounts for the discreteness of the particles and is assumed to obey Poisso- 
nian statistics with 



;e(x)) = o 

(£(x)t(x + A)) = {N(x))8(£). (2.14) 

Note, that N(x) is a measured quantity and not a quantum mechanical operator or an 
ensemble average. This implies that both N(x) and £(x) may differ between individual runs 
of the experiment. Since the individual noisy contribution £(x) is unpredictable, all we can 
hope to achieve, is a complete characterization of the smooth patterns N(x). 

It was indicated above, that an analysis of correlation functions will be required to make 
substantial statements about the behaviour of single runs. The lowest order correlation 
function is given by 

corrA^(x, A) = (n(x)N(x + A)), (2.15) 

where the angular brackets denote an average over several runs of the experiment. We can 
calculate this correlation by identifying the measured density distribution N(x) with the 
corresponding quantum mechanical operator ^^(x)^f(x). This yields 

corrA^(x, A) = (&(x)V(x)&(x + A)*(ac + A)). (2.16) 



Since we aim to characterize the smooth patterns N(x), we are interested in correlations 
of N(x), defined analogously to Eq. (|2.15| ), instead of those of the measured density N(x). 
Therefore, we rewrite Eq. ( |2.15|) using Eq. ( |2.13|) and Eq. (p. 14 ) as 
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corriV(x, A) = corriV(x, A) + (£(x)£(x + A)) (2.17) 
= corriV(x, A) + (n(x))S{A). 



Normal ordering of the field operators in Eq. ( [2.16 ) generates a term that is identical with 
the delta correlated noise in Eq. ( 2.17Q . The correlation function for the smooth pattern is 
thus given by 

corriV(x, A) = (^(x)^(x + A)^(x + A)#(x)Y (2.18) 



We will use this equation in the following to determine the distribution of all smooth patterns 
N(x) that can be seen in single runs of the experiment. 

Since the condensates are described as macroscopically occupied momentum modes, we 
have to link the above description in position space with one in momentum space. The 
momentum decomposition of the quantum mechanical field operator along the x-axis on the 
screen is 

*(x) = -jL= y e lkx a k (2.19) 

while any possible smooth interference pattern can be described by the Fourier decomposi- 
tion 

N(x) = -^J2e ik *h k . (2.20) 

k 

Note that h k is a complex random variable with yet unknown statistics. It satisfies the 
relation n k = n*_ k to ensure that N(x) is real. In the following we will derive the statistics 
of the Fourier components n k by a comparison of the higher moments of h k with higher 
moments of the quantum mechanical operators a k . 

Identification of the average over the measured screen patterns (iV(x)) with the quantum 

mechanical prediction (ty(x)^(x) 



rkrs E z lkx = T^- 3 E <? kx E (2-21) 

k k k' 

yields 

(^> = E(4- fc «fc') (2-22) 



due to the uniqueness of the Fourier decomposition. Since the two momentum modes are 
assumed to be independent we get 

(n k ) = 4,o (iVi + N 2 ) . (2.23) 

Note, that the ensemble averages (a k ) vanish. In a similar way the correlation of the mea- 
sured patterns corriV 
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corrA>(:r, A)V 2/3 = ^ e i(fc+fc,) V fcA (h k n k ,) (2.24) 

kk' 

can be identified with its quantum mechanical analogon 

corriV(a;, A)V 2/3 = (a\a\aiai) + {a\a\a 2 a 2 \ (2.25) 
+2A^iA^ 2 (l + cos2A;oA). 



Eq. ( p.25 ) states that the probability for observing an atom at a position x + A after 



having observed a previous atom at position x is a periodic function of 2A; A. Already this 
observation implies that the density distribution of a single run must show some spatial 
variation. Eq. ( [2.24j ) and Eq. (|2.25| ) determine the only nonvanishing coefficients {hkhk') 



(h 2 ) = (aXala^) + (44a 2 a 2 ) + 2AqiV 2 , (2.26) 

2^ 



\n ±2ko \ 2 ) = N,N 2 . (2.27) 
In the case of small number fluctuations, as assumed in Eq. (|2.8|) , we get 

(njj) = (JVx + N 2 ) 2 = (n f . (2.28) 



This equality implies that n$ does not fluctuate. Rather than being a random variable, its 
value is invariably given by 

n = N x + N 2 . (2.29) 



Eq. (|2.27|) and the vanishing expectation value of n 2 k Eq. ( |2.23| ) show that at least the 



phase of the Fourier coefficient n 2 k is subject to fluctuations. A similar evaluation of the 
fourth order correlation function yields 



h 2 kX) = (N 1 N 2 ) 2 . (2.30) 

We skipped the explicit calculation here. One sees by a comparison with Eq. (|2.27|) that the 
modulus of h 2 k does not fluctuate either. It is fixed to the value 



|n 2fco | = V^2- (2.31) 

Since the classical correlation functions are of even order in the atomic field operators, any 
phase factors cancel and no constraint for the phase of h 2 k exists. 

Thus, Eq. ( p. 29 ) and Eq. (j2.31|) provide us with the complete obtainable information 



about the Fourier coefficients of the interference patterns (cf. Eq. (|2.20|) ) that can be mea- 
sured in single runs of the experiment. Since no constraint for the phase of h 2 k exists, its 
value if is an equally distributed random variable. The result of this derivation is, that the 
distribution of all possible recorded smooth interference patterns N(x), rather than the full 
noisy patterns N(x), is given by Eq. fl2.11|) . 



A surprising consequence is that with Eq. ( p.ll| ) holds (up to order 1/iVj) no matter 
whether the density matrices of the individual modes are given by Eq. fl2.6|) or Eq. ( p.7[) . 
This shows that a Fock state will yield in the limit of large particle numbers virtually the 
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same experimental results as an incoherent superposition of coherent states (see also Ref. 
H). The neglected terms of order 1/iVj mainly account for the existence of some correlations 
in the earlier eliminated noise £(x) in order to preserve the total particle number precisely. 

We close this section with a few remarks about the range of validity of the above result 
under more realistic conditions. We have assumed in our calculations that the fluctuations 
in the flux are negligible during the detection time r^. Inclusion of these fluctuations will 
clarify the range of validity of the above considerations. When taking into account that the 
flux of particles may vary during the formation of the interference pattern, the correlation 
function Eq. Q2.15| ) changes to 

corriV(x, t; A, r) = (n(x, t)N(x + A, t + r)). (2.32) 



Eq. ( 2.25Q can then be rewritten in terms of the time correlation functions that characterize 



the dynamics of the two independent modes (i = 1,2) 
corriV(A, t)V 2/3 = N*g{ 2 \r) + N 2 g { 2 2) (r) 



+2N X N 2 (l + cos2k ARe {g^ (^g^* (r)}) (2.33) 



with 



g?\r) = ( Qt(t)a(t + r) / ^ Q 

y/(at{t)a{t)) (at{t + T)a{t + T)) 

™ (aUt)aUt + r)a(t + r)a(t)) 

9?\r) = }„\ , \, J z { v ™ 1. (2.34) 

yi w (at(t)a(i)) (at(t + r)a(t + r)) v ; 

It is evident that the decay of the single interference pattern during a longer exposure time 
is intimately connected to the decay of the amplitude correlation function g^{r). Such a 
decay will happen in any finite physical system. g^(r) will decay on a similar time scale 
and we get 

lim corriV = (Aq + N 2 f . (2.35) 

T— >00 

This implies that even a single run is appropriately described by the quantum mechanical en- 
semble average (ty^(x)ty(x)\ as we are used from ordinary atomic interference experiments. 
However, in a system with a strong quantum degeneracy like a laser or a Bose condensate 
the decay of g^(r) may be so slow compared to the flux of atoms so that intermediate time 
scales are relevant for experimental purposes. The concept of broken symmetry applies only 
for these time scales. 



III. INTERFERENCE OF TWO LOCALIZED CONDENSATES 

In the last section we approximated the condensates by two plane waves. The general 
result, however, is equally valid for more realistic setups. Here we show by a numerical 
simulation that the setup of || indeed implements such an interference experiment, although 
at present no interference has been observed due to an insufficient experimental resolution. 
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The experiment, reported in ||, can be summarized as follows. Two independent Bose 
condensates of 23 Na atoms are stored in a magnetic trap, separated by a laser beam. These 
constraints can be modelled by two identical harmonic potentials with angular frequencies 
v = 2tt x 745s -1 , 2tt x 235s~ 1 , 2tc x 410s~ 1 in the x,y, and z directions, respectively. These 
potentials are separated along the x-axis by xo = 100/im. The total number of atoms is 
assumed to be Nq = 1.5 x 10 5 , or 0.75 x 10 5 atoms in each condensate. The scattering 
length, characterizing the atomic interaction, amounts to a = 4.9nm. After switching off 
the constraints, the atomic clouds expand and eventually overlap, both due to the quantum 
mechanical spreading of their wave packets and to the atomic repulsion. The final spatial 
distribution of the atoms is measured after 6ms or equivalently 28i , ~ l . 

In our calculation we simulate the expansion of these atomic clouds. We assume that the 
two condensates initially are at zero temperature, with both condensates equally populated. 
The investigation of atomic interference at nonvanishing temperatures is put forward to the 
next section. There, it is shown that interference can be observed in a different setup also 
at small but finite temperatures. 

The sodium atoms are subject to a repulsive interaction. The initial wave function of a 
single condensate is therefore described by a solution of the Gross-Pitaevskii equation || 

^(r,t)= |-^ + ^ + W(r,*)| 2 | ^(f, t) (3.1) 



with the chemical potential /x and 



U = N^. (3.2) 



m 



Note that the macroscopic wave function if)(f,t) is normalized to one in this notation. The 
chemical potential 

U ( N \ 2/5 

» = - m T {te*wlz*aW) (3 ' 3) 

is expressed in terms of the natural length units of the harmonic potential V(f) along the 
coordinate axes. In the case of the x-direction, for example, this length unit is 



Ax = \ = 0.78/mi. (3.4) 
V fnv x 

Since the total energy is dominated by the particle interactions and the external harmonic 
potential, it is a good first approximation |T6j to neglect the kinetic energy in Eq. fl3.1| ). This 



yields the approximate expression 

ij(f) « J(n-V(ij)/U (3.5) 

for the initial wave function, provided that the expression under the square root is positive. 
Otherwise, the wave function vanishes abruptly, for example on the x-axis at a distance 
of 3.24/iin from the origin. The physical wave function, however, vanishes smoothly on a 
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small scale that is called the healing length |16|. This smooth rather than abrupt vanishing 



accounts for the existence of a limited kinetic energy. 

We therefore conclude that the initial wave function of a single condensate is largely 
of a parabolic shape and about 4 times larger than the corresponding noninteracting wave 
function. Note, that this significant deviation from the pure harmonic oscillator ground 
state is due to a balance between the potential energies of the atom interactions and the 
external harmonic potential. It should not be confused with the depletion of the condensate 
that is found in an interacting homogenous gas M: 



iV-iVn 8 na 3 , , 

(3.6) 



N 



7T 



Inserting the peak densities of the current experiment into this formula indicates a depletion 
of less than 1% which is negligible for our purposes. 

According to the results of the last section it is justified to describe the two independent 
condensates by one combined macroscopic wave function, for characterizing a single run of 
the experiment. This wave function consists of two spatially separated parts, corresponding 
to the two condensates, that are joined with a relative phase if between them. Of course, 
this relative phase varies between different runs of the experiment. Since variations in the 
relative phase lead only to trivial shifts of the interference pattern we restrict ourselves to 
the case ip = 0. 

The trap is switched off at a given time and the atomic clouds expand. We describe this 
expansion assuming the validity of the nonlinear Schrddinger equation [16| 



ih^(r,t)= l-^ + U\ij(r,t)\ 2 ]ij(r,t). (3.7) 



This is also referred to as the time- dependent Gross-Pitaevskii equation (cf. Eq. (|3.1|) ). 
Its formal structure resembles a mean-field equation when replacing the original two-body 
potential by a pseudopotential. It can be derived for cold and dilute gases by applying the 
ladder approximation to the two-body Green's function |Q. Even within this approximation 
Eq. ( |3.7|) does not account for the full dynamics, but it should be a good approximation for 
time scales, where collisional rates can be neglected.0 

Due to the nonlinearity in Eq. ( |3.7|) , it is incorrect to separate the wave function with re- 
spect to the spatial coordinates. Nevertheless, a separation of the spatial degrees of freedom 
has to be applied as an approximation, since a full three-dimensional numerical treatment 
of the problem is far beyond reach. Such an approximation goes back to |17] and has found 



widespread use in chemical physics where it is called time-dependent self-consistent field 
(TDSCF) method. Its validity has been studied for example in ||18|| . This approximation 



reduces the propagation of a three-dimensional wave function to that of three coupled one- 
dimensional wave functions. The error of this simplification depends on the deviation of the 
wave function from a Gaussian distribution. 



Note, that the atomic clouds are not in thermal equilibrium during the described expansion. The 
earlier assumption of zero temperature applies only in the presence of the trapping potentials. 



12 



The numerical propagation of the wave function is performed by a split operator tech- 



nique that is accurate and easy to implement. The initial wave function is generated 



by an imaginary time propagation of Eq. ( |3.7| ), which yields the full solution of Eq. ( |3.1| ) in 
contrast to the approximate wave function Eq. ( |3.5|) . 

We plot in Fig. [I] the time evolution of the x distribution, as it could have be seen in 
the experiment. We have averaged over the relative phase of the two condensates since the 
fringes are not resolvable on the plotted scale. The two wave functions overlap after about 
3ms. Even though there is a considerable overlap of the two condensates, it is not large 
enough to allow for an interpretation of the x distribution in terms of a velocity distribution 
of the initial single condensates. 

In Fig. ^| we show both the initial and final x distribution including the interference 
fringes as they are predicted by the nonlinear Schrodinger equation. The periodic length 
of the fringes is about 1/im and thus below the pixel resolution of 12.4/im used in the 
experiment of [[|. The period of the fringes increases slowly with time. However, its size is 
of the same order of magnitude as the healing length throughout the whole investigated time 
interval. The period of the fringes could be enlarged by widening the trap potentials before 
dropping the atoms. This would reduce both the initial kinetic and interaction energy. 

Fig. |3| shows the calculated x-z distribution after 6ms. Its square like shape is an ar- 
tifact of the TDSCF approximation. A full three-dimensional calculation would yield a 
slightly larger extension along the axes. For comparison, we plot in Fig. |] the corresponding 
distribution, calculated from the experimental data of 0. The original data describe the 
spatial dependence of the light absorption. Assuming an exponential absorption law we 
have transformed the light absorption into a particle distribution and normalized it for a 
better comparison with Fig. [| We conclude that the experimental data of || are in good 
agreement with a description using the nonlinear Schrodinger equation. While the existence 
of shoulders and their positions are in especially good agreement, also the total size of the 
atomic cloud is fairly well reproduced. A detailed quantitative analysis, however, requires a 
better accuracy of the experimental data. 



IV. INTERFERENCE PROPERTIES OF A SINGLE BOSE CONDENSATE 

In this Section we analyze the effects of finite temperatures in interference experiments 
with a Bose gas. We consider a sample of N bosons confined in a three-dimensional harmonic 
trap, not necessarily isotropic. After reaching thermal equilibrium at a given temperature 
T, the gas is released and the sample falls through two pinholes, separated by a distance 2d, 
that are located on a screen Si at a distance L from the trap center (Fig. The atoms 
are recorded on a second screen S2, at a distance D from the trap center. We wish to stress 
the fact that in this experiment the wavefunctions coming from both pinholes have fixed 
relative phases since they represent filtered parts from the same large wavefunction that was 
stored in the trap. Therefore, different realizations of the experiment will show the same 
interference pattern. This is in contrast with the experiments analyzed in Sec. 0, in which 
we considered two independent condensates. We will not perform the correlation analysis of 
Sec. f| here since it does not yield additional information in this case. The main purpose 
of this Section is to show how the interference pattern depends on the temperature of the 
sample, and to what extent it reflects the phase transition at the critical temperature. We 
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will consider the case of an ideal gas, i.e. we will not take into account the role of atomic 
collisions. The problem of a (weakly) interacting Bose gas at finite temperatures deserves a 
separate analysis. 



A. Evolution through the pinholes 

We are interested in the mean number of particles dN per unit time and unit area 
deposited at any point f = (x, y, D) on the screen 5*2 at a given time t. This quantity is 
given by the expectation value of the z component of the probability current, 

where M is the atomic mass, and \l/(r, t) is the field operator describing the sample. For the 
case of non-interacting particles, this operator can always be written as 

*(f,t) = y E,a i pl> i t(r,t), (4.2) 

n 

where the wavefunction ipft(r, t) satisfies the (single particle) Schrodinger equation describing 
the evolution of a single particle. Here, ^(f, 0) is the eigenfunction of the three-dimensional 
harmonic oscillator with quantum numbers n = (n x , n y ,n z ) (n — 0, 1, . . .). As usual, are 
annihilation opeartors of particles in the state n. 

It is worth mentioning that, if the wavefunction is expanded in terms of momentum 
eigenstates as in the previous Section, I(f,t) can be fully characterized by the density of 
particles. Here, however, we find it more convenient to use the eigenstates of the harmonic 
oscillator. Therefore the probability current has to be evaluated explicitly for describing the 
interference pattern. 

Substituting (|4.2|) into (|4.1|) , we obtain 



I(r,t) = Y^N n I H (r,t) (4.3) 

n 

implying that the incoherently populated trap eigenf unctions propagate independently after 
releasing the trap. Their occupation numbers = (atag) are given by the Bose-Einstein 
distribution 

\ p -f3h(n x v x +n y u y +n z v z ) 

= (4.4) 

n \ \q— fih(n x Ux+n y Uy+n z u z ) \ ' 

where (3 = 1/kbT, A is related to the number of particles in the trap N, and v x ,y,z the trap 
frequencies along the three axes (see Appendix B). We have also defined 



T t- \ ^ T 

Ifi(r, t) = — lm 



d 

$ft{r,t) —ij;fi(r,t) 



(4.5) 



In view of this formula, we only need to derive the expression for ^(r*, £); that is, we simply 
have to find the evolution through the pinholes of a single particle that is initially in the 
harmonic oscillator eigenstate with quantum numbers n. 
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Apart from neglecting the role of collisions, in order to calculate the evolution of if}a(f, 0) 
we will make other simplifying assumptions. First, we will assume that L is small enough 
such that during the time that the sample needs to reach the pinholes, this wavefunction 
basically does not change. Second, we will assume that the size of the pinholes is small 
compared to the typical distances over which this wavefunction varies. Finally, we will take 
D to be large enough in such a way that the parts of the wavefunction coming from each 
pinhole overlaps with each other. According to these assumptions, we can approximate 



1pn(r,0) ~ klp nx (x)[5(x - Xi) + 6(x - X 2 )]lpn v (y)S(y - yo)lpn z ( z ), 



(4.6) 



where (xi, yo) and (x 2 , yo) are the coordinates of the pinholes in the first screen Si, ip n (x) is 
the eigenstate of the one-dimensional harmonic oscillator in position representation, and k 
a normalization constant. 

After a time t > 0, the wavefunction evolves to 



ipn(r,t) ~ k[G (x,x u t)i/) nx (xi) + G (x,X2,t)il) nx (x 2 )}Go(y,y ,t)il) ny (y )tlj nz (z,t), 
where 



Here, 



G (x,x',t) = e~ in/i 
G g (z,z',t) = e-^ A 



dz'G g (z,z' ',t)ijj nz (z'). 



(4.7) 



(4.8) 



M 



X l 2 .M(x-x>) 2 



27iht 

r M ] V2 e i^- t [(z-z') 2 + 9 t 2 (z + z')-g 2 t^ 

2nht\ 



(4.9a) 
(4.9b) 



are the free propagator and the propagator under a constant force, respectively. The sub- 
stitution of Q4.7| ) into fl4.5[ ) allows us to write 



where 



l n (r,t) = k*Il(x,t)P ny (y,t)I* z (z,t) 



r nS x ^) = |Go(x,Xi,t)^(a;i) + G (x,x 2l t)^ nx (x 2 )\ 2 , 

^(y>*) = \ G o(y,yi,t)ipn y (yo)\ 2 , 



h 



:Im 



d 



(4.10) 



(4.11a) 
(4.11b) 

(4.11c) 



Let us now derive simple expressions for these quantities. First, using (|4.9a| ) we find 

n (y,t) = (M/2nht)\^ ny (y )\ 2 (4.12) 



for all y [20]. On the other hand, as it is shown in the Appendix A, 

r n (z,t) = F{t)\i, nz {z - gt 2 /2,0)| 2 , 



(4.13) 
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where 



F(t) 



gt 

i + {v z ty 



,2 r- 



1 + 



(4.14) 



and ip nz is the eigenstate of the one- dimensional harmonic oscillator with frequency u z /[l + 
(v z t) 2 ]. Finally, 

Il(x,t) = (M/27iM){\ij nx (x 1 )\ 2 + |^(x 2 )| 2 + 2^(x 1 )W(z 2 ) cos[2nx/x f (t) + 0(t)]} 

(4.15) 

where we have taken into account that ip nx (x) is real [pOfl . Here, 



27(f) 



7T?l£ 



2'KV T t 



''Ox 



(4.16) 



with aox = [?i/(2M^)] 1 / 2 being the size of the ground state wavefunction, and 4>(t) 
M{x 2 -x 2 )/{2M). 



B. Interference and visibility 

Let us now specialize the above derived expressions to a simple case. We take yo = and 
x\ = d = —X2 (i.e. the pinholes are symmetrically situated along the x axis). For the sake 
of simplicity we choose the time r = {2D / ' g) 1 ^ 2 , which is the time required for the center of 
the trap to reach by free fall the screen 5*2. In this case, we obtain 

/*>,r) cx |</Urf)| 2 {l + (-l) n *cos[2nx/x f (T)}}, (4.17a) 
^>,r)oc|^(0)| 2 , (4.17b) 
/:>,r)och/U0)| 2 , (4.17c) 

Note that the term (— l) nx appearing in the first of these expresions comes from the fact 
that ip nx {—d) = (—l) nx ip nx (d). Thus, depending on n x being even or odd, the interference 
pattern has a maximum or a minimum at x = 0. It is obvious that the incoherent addition of 
all the contributions for different quantum numbers n x [see ( [Op ] will decrease the visibility. 
Therefore, the origin of the temperature dependence of the interference fringes is precisely 
the term (— l) Ux accompanying the cosine. For sufficiently low temperature, the Bose- 
Einstein distribution ( |4.4| ) is peaked at n x = 0, and therefore fringes along the y axis will 
be clearly displayed on the screen S*2. The separation of the fringes is Xf(r) given in ( |4.16|) . 
Note that, in principle, this separation changes with time. In practice, this variation will be 
negligible as long as the time required by the sample to cross the screen 5t is much smaller 
than r. More specifically, we can neglect this time dependence for 

y = vMl + M 2 ]|§ « 1. (4.18) 

This condition can always be fulfilled for a sufficiently long distance D (note that r scales 
as v^D, and therefore 5t/t scales as 1/y/D). 
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Let us now analyze the temperature dependence of the visibility 



rmax rmin 

V=- — . (4.19) 

Jmax _|_ Jmm v ' 



where I is given in ( ]4.3|) . The maxima and minima of the expression ( |4.17a| ) can be easily 
calculated, taking into account that for fixed n VtZ , is a decreasing function of n x . Thus, 
the maxima (minima) in the fringes correspond to x™ ax = kxf(r) [x™ in = (k + l/2)xf(r)], 
with k = 0, ±1, . . .. Using this result we find that 

mr^^^i^MiV^onVn^o)! 2 , (4.20a) 

n 

jmax _ jmm ^ ^(-l)^^|^(d)i 2 |^(0)| 2 |^(0)i 2 . (4.20b) 

n 

In Fig. |]we have plotted the visibility (solid-lines) as a function of the scaled temperature 
kbT / (his) for N = 10000 particles in an isotropic harmonic oscillator and different distances 
2d between the two pinholes. We have also plotted (dotted-line) the proportion of particles in 
the condensate N (for the numerical method used to calculate these results see Appendix 
B). Note first that for a given temperature, the visibility decreases as d increases. The 
reason for that is as follows. As mentioned above, the eigenstates of the harmonic oscillator 
with even quantum numbers n x give in phase contributions to the fringe pattern, whereas 
the ones with odd n x give out of phase contributions, i.e. tend to decrease the visibility. 
The effect of the terms with odd n x tends to zero for d — > 0, since ip nx (d) tends also to zero 
(note that this wavefunction is antisymmetric). Thus, the pinholes select the particles whose 
wavefunctions have even n x . As d increases, up to the order of the ground state wavefunction 
ao x , the visibility decreases since the contributions of the odd wavefunctions becomes more 
important. For d 3> a 0x the visibility remains small even below the critical point, since only 
the wavefunctions with high n x contribute. On the other hand, the transition point of BEC 
is reflected in the visibility, which shows that double-slit experiments can be used to observe 
the phase transition. Note the dramatic changes at the transition point for moderate values 
of d. It is somehow surprising that even for N — No > Nq the visibility already becomes very 
large, since one would expect that the particles that are out of the condensate would produce 
some noise that would dominate over the effect of the condensed particles. However, this 
is not the case. The N — N particles out of the condensate are distributed among several 
states with different n x (this distribution is similar to a Boltzmann distribution and therefore 
has a long tail). As n x increases, their contribution to the intensity becomes smaller and 
smaller given that \ipn x (d)\ 2 ~* 0. Thus, for d ^ ao x the change of visibility at the transition 
point becomes very pronounced. 



V. CONCLUSION 

We have shown that two independent Bose condensates, when overlapping on a screen, 
can exhibit an interference pattern in a single run of the experiment. A rigorous analysis in 
terms of correlation functions confirmed that the description of an individual Bose conden- 
sate by a coherent state, i.e. a macroscopic wave function, gives the proper characterization 
of single runs of this interference experiment. This applies although atomic coherent states 
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violate fundamental conservation laws and holds in the limit of large occupation numbers 
even for a single Fock state. The macroscopic wave function concept resembles the way 
Quantum Electrodynamics converges to Classical Electrodynamics for large photon num- 
bers. However, one fundamental difference between atomic and photonic coherent states 
should be mentioned. It is, at least in principle, possible to measure the absolute phase 
of a photon coherent state, since the electric field amplitude is a true observable. As a 
consequence, the photon number is no conserved quantity. In contrast, the atom number is 
strictly conserved. In return, the macroscopic wave function itself is not an observable, thus 
inhibiting the measurement of its absolute phase. Consequently, the investigated interfer- 
ence experiment yielded only information about the relative phase between two macroscopic 
wave functions. 

It is another aim of the paper to describe simple experiments that allow for a comparison 
between experimental and theoretical features of a Bose condensate. The considerations 
above, concerning the interference of independent Bose condensates allowed the numerical 
simulation of a recent experiment where two independent Bose condensates are released 
from a magnetic trap. They expand and eventually overlap. In our calculations we accounted 
for particle interactions but assumed an initial distribution with zero temperature. The 
calculations predict interference fringes with a period below the resolution of the current 
experimental data. However, the size of the final atomic cloud is in good agreement with 
the experimental result. A further quantitative analysis is limited by the present accuracy 
of the experiment. By improved detection techniques it should be possible to study the 
existence of the interference fringes. This could yield substantial and new information about 
the condensate and the range of validity of the nonlinear Schrodinger equation. 

As a possible future experiment we proposed the study of the interference properties of a 
Bose condensed system in a conventional Young interference experiment. We assumed that 
the atoms are released from a harmonic trap and fall freely through two pinholes. Here we 
neglected particle interactions since the above calculations have shown that the visibility of 
the interference fringes is not severely affected by these interactions. Instead we studied the 
dependence of the visibility on the temperature. The transition to a Bose condensed phase 
is reflected by a sharp rise in the visibility of the fringes. Under favourable conditions the 
visibility reaches values close to one immediately below the critical temperature. In addition 
a variation of the distance of the the pinholes shows that the coherence length of a finite 
Bose condensed system is comparable with the size of the groundstate wavefunction. 
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APPENDIX A: EVOLUTION IN THE Z DIRECTION 



In this appendix we derive some of the formulas used in Section III. First, starting from 

(Al) 



(|4.8|) we write 

/OO 
dz'H n {az')e- mz ' 2+tC{z ' t)z ' 
-oo 

Here, we have used 

<4, n {z^) = K n H n {az)e- a2z2 l\ 



(A2) 



where a = (Mv z /h) l l 2 , K n = [a/ '(2 n n\^7r)] 1 / 2 , and H n is the n-th Hermite polynomial. We 
have also defined 



A(z,t) = 




/4 / M JM/{2M){z 2 +gt 2 z-gH*/l2) 

V 27tM 


C(z,t) = 


M 
M 


r 1 2 

-gt - z 
2 


j 


R(t) = 


a 2 
2 


. M 
~ l 2M ~ 


M , 

-—-\v z t — i). 



(A3a) 
(A3b) 
(A3c) 



Performing the change of variables J R{t)z = x and transforming the path for the integration 
in the complex plane, we arrive at 



A (z,t) 
R(t) 



K n e~ B ^ 2 / dz'H, 



az 



R(t) 



[x-iB(z,t)] 2 



(A4) 



where B(z,t) = C(z, t)/[2yR(t)]. Performing the integration, and after some lenghty alge- 
bra we obtain 



A (z,t) 
R(t) 



Kne -a(trz( z ,tf/2[i+i/(v z t)] l- R (t)*/R(t)} n/2 H n [&(t)Z(z, t)}, (A5) 



where a = a/^J 1 + u 2 t 2 and Z(z,t) = z — gt 2 /2. Now, the expression (|4.13|) can be easily 
derived starting from ( [4.1 lc| ) and using (|A"5|). To do that, one has to note that when 
performing the derivatives of if} n (z, t), the contribution given by the derivative of the Hermite 
polynomial vanishes when taking the imaginary part. We obtain 

P nz {z,t) = F{t)\^ nz {z,t)\\ (A6) 
where F(t) is given in ( [4.141) . Taking the modulus square of ( |A5f) we obtain (|4.11c|) . 
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APPENDIX B: FORMULAS FOR THE NUMERICAL EVALUATION OF THE 

VISIBILITY 



In this appendix we give some of the formulas used to evaluate numerically the expression 
of the visibility for a Bose-Einstein distribution. The total number of particles can be written 

as 



n k=l 



iVn 



(Bl) 



i=x,y,z 



Given a fixed number of particles N and a temperature T, we have first determined the 
value of N Q by using a bisection method varying N until Eq. ( |Blp is verified. Thus, we can 
calculate the whole distribution through Eq. (|4.4j) since A = N /(l + N ). 

On the other hand, in order to calculate the expressions ( |4.2U| ) we need to determine 
sums of the form 



Sl = J2e- PhvHn+1/2) \Mx)\ 2 , 



n=0 
oo 



S 2 = ^(-l)" e - /37il ' fc(n+1/2) |^ n (x) 



(B2a) 
(B2b) 



n=0 



This can be easily performed by noting that these expressions are related to the propagator 
for the harmonic oscillator with an imaginary time. In particular, 



51 = G(x,x,t = —i(3k) 

5 2 = G[x, x, t = tt/ {hv) — i/3k], 



where 



G(x,x,t) = ^A/{2ix)e Ax2[l - cos{vt) \ 
with A = Mi//[ihsm(i/t)]. Using these expresions, we find 



rmax , rmin 



+ r 



x N - k 



rmax rmm 




,« 2 d 2 [l-cosh(/3hu x k)] I smh(/3hv x k) 



k=l 



Nq + 1. 



i=x,y,z 

-a 2 d 2 [l+cosh([3hv x k)]/ sinh(f3hu x k) 

i=x,y,z 



Yl ^smh(Phuik), 



(B3a) 
(B3b) 



(B4) 



n v smh (^ fc )' ( B5a ) 



(B5b) 
(B5c) 



Note that one could use directly the formulas (f4.20|) , which involve 3 nested sums. However, 
Eq. (|B5|) involves a single sum, which saves a lot of computer time. 
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FIGURES 



FIG. 1. Calculated time evolution of the x 
distribution. The interference fringes are re- 
moved by averaging over <p. The contour lines 
correspond to an increment of 0.001. 

FIG. 2. Initial and final x distribution. The 
final distribution is plotted both with a rela- 
tive phase if = and averaged over all possible 
phases. 

FIG. 3. Calculated final x-z distribution 
with removed interference fringes. 

FIG. 4. Experimental atomic density, calcu- 
lated from the light absorption. A renormaliza- 
tion was performed to facilitate comparison with 
Fig. |3|. The original data correspond to Fig. 2c of 
||. Their use is with kind permission of W. Ket- 
terle and coworkers. 

FIG. 5. Setup of the experiment considered 
in Section III. An ideal Bose gas in thermal equi- 
librium at temperature T is dropped through 
two pinholes. The interference is recorded on 
a screen. 

FIG. 6. Visibility V plotted as a function of 
the scaled temperature for (curves from top to 
bottom) d = 0.5,1.56,2.61,3.67,4.72,5.78, and 
6.83[h/(Mv)] (the inset shows a detail). In dot- 
ted line we have plotted the proportion of parti- 
cles in the condensate no = Nq/N. 
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